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ABSTRACT 

Fermi observations of GeV emission from GRBs have suggested that the 
Lorentz factor of some GRBs is around a thousand or even higher. At the same 
time the same Fermi observations have shown an extended GeV emission indi- 
cating that this higher energy emission might be a part of the afterglow and it 
does not come from the same region as the lower energy prompt emission. If 
this interpretation is correct then we may have to reconsider the opacity limits 
on the Lorentz factor which are based on a one-zone model. We describe here a 
two-zone model in which the GeV photons are emitted in a larger radius than 
the MeV photons and we calculate the optical depth for pair creation of a GeV 
photon passing the lower energy photons shell. We find that, as expected, the 
new two-zone limits on the Lorentz factor are significantly lower. When applied 
to Fermi bursts the corresponding limits are lower by a factor of five compared 
to the one-zone model. It is possible that both the MeV and GeV regions have 
relatively modest Lorentz factors (~ 200 — 400) , which is significantly softer then 
one zone limit . 

Subject headings: gamma rays: bursts — radiation mechanisms: nonthermal - 
relativity 

1. Introduction 



Highly relativistic motion, essential to overcome the Compactness problem ( IRuderman 

19751 ). is a basic ingredient of all GRB models. The value of the bulk Lorentz factor, T, 
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of the relativistic outflow is of outmost interest. It is essential for understanding the na- 
ture of the inner engine, the outflow and its acceleration and collimation mechanisms, the 
conditions at the emitting regions and the radiation mechanism. So far the most robust 
method to estimate Y was using the Compactn ess. The high energy photons s et an upper 



limit on the optical depth for pair product ion (IFenimore. Epstein fc Holll993t iPiranl 11995 



Woods fc Loeblll995l ; iLithwick fc Sarill200ll ). The observations of GeV photons from several 
bursts enabled the Fe rmi team to set very hig h (r > 1000) lower limits on T for those GRBs 
JAbdo et al.ll2009al lbl: kckermann et aliboioh. While other methods to estimate or limit T 
depend on various assumptions (see e.g. IZou &: Piranl (120101 )). the compactness limit seems 
to be independent of any model assumptions and hence it is considered to be the most robust 
one. 

However, recent Fermi observations have also shown that the onset of higher energy 
(~ 100 MeV to a few GeV - denoted hereafter GeV) emission lags after the onset of the lower 
energy ( ^100 keV to a few MeV - denoted hereafter MeV) prompt emission (jAbdo et al. 
2009al lbl; lAckermann et al" 2010 ; Ghisellini et al. 2010 ). Fermi confirm ed earlier EGRET 



results that the GeV emi s sion also lasts longer than the MeV emission (IHurley et al.l 11994 



Fishman fc Meeganlll995l ; iGonzalez et al.l 120031 ) . These two facts suggest the possibility of 
a different origin for the MeV and the GeV emission. This would be the case, for example, 
if the GeV emission arises from an external shock afterglow (Kumar & Barniol Duran 2009, 
2010; G hisellini et al. 2010; P i ran & Nakar 2010; Gao et al. 2009), f rom multi-zone internal 
shocks JXue. Fan fe Weill2008l : koi et alJliooa Izhao. Li fe Bafeoioh or if the MeV emission 
is the quasi-thermal radiation of the baryonic outflow while the GeV emission is mainly from 
the subsequent internal shocks^. 

The compactness limits on the Lorentz factor are based, however, on an implicit as- 
sumption that the MeV and the GeV photon arise from the same region. We show here that 
the relaxation of th is a ssumption reduces significa ntly the estimated lower limit on Y (see 
also lAoi et al.l |2009| ; iLil I201CH : IZhao. Li fc Bail 120101 ) . The existence of two regimes leads to 



a rich variety of possibilities. We consider in the following the most natural configuration, 
which is also consistent with the temporal delay of the GeV emission, the MeV emission 
is produced at lower radii (say via internal shocks) and the GeV emission is produced at a 
larger radius (say via external shocks) - see Fig. [TJ Other geometrical options that we don't 
consider are that the GeV emission is emitted at a lower radius than the MeV emission or 
that the MeV and the GeV emission are produced at different angular regions. Our estimates 



lr The synchrotron radiation of electrons accelerated by such internal shocks may peak at eV energies while 
the inverse Compton radiation can give rise to a significant GeV component. In such a scenario, the GeV 
emission and energetic ultraviolet/optical flare are tightly correlated. 
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don't depend on the origin of the emission (e.g. internal or external shock) but just on the 
overall geometry of the system. We calculate the optical depth of a GeV photon passing 
through the MeV photons shell and we obtain new compactness limits on T. 



2. The Model 

Consider two emitting regions denoted by M, for MeV and G, for GeV respectively (see 
Fig. [1]). The MeV (GeV) emission region has a radius Rm {Rg) (for simplicity we consider 
emission from thin shells) with R M <C Rg, a Lorentz factor T M (T G ) and an angular width 
@m,j {9g,j)- We assume that the MeV and GeV jets are aligned and both are pointing toward 
us. In this configuration a GeV photon passes through a shell of MeV photons on its way to 
the observer (see Fig. [2]). Our goal is to estimate the optical depth for pair production of a 
GeV and an MeV photons shell. The width of the MeV photons shell ii^. 

A M ~cT m /(l + z), (1) 

where T 90 is the observed duration of the MeV pulse and z is the redshift of the burst. A 
GeV photon, emitted along the x axis at R G and 9 G is immersed in MeV photons until it 
leaves the MeV shell at R m ^( R G,0 G ) (see Fig. [2]): 

Rg 

R^Rg, 9 G ) = 1 _ 2a > Ml + z) 02 oRQ y ( 2 ) 

for 9% > 2cT 90 /[(l + z)R G ], while R m ^(R G , 9 G ) -+ oo for 9% < 2cT 90 /[(l + z)R G }. 

Assume, for simplicity, that the emitted flux of the MeV photons is constant over time 
then the number density of the MeV photons is: 

U0 > mMR) " (l + zyE p AnRV (3) 

where L M is the isotropic equivalent MeV luminosity, E P is the observed peak energy and 
R G < R < Rmax (R G , 9q) is the radius. The spectrum of the MeV photons is described by 
the Band function: 

' '§?)"" E < E - 

n Me v(E M ) ^ n 0MeV (R) { ) : ' Pm (4) 

§^-J E p < E < i?M,max, 



IZhao. Li fe Bail (|2010h consider erroneously an MeV shell width of R/2T (see their Eq. (10)). 
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GeV region 

Fig. 1. — A schematic diagram of the two-region scenario. GeV photons are produced at a 
larger radius. The MeV shell engulf this radius when the GeV photons are produced. At 
Rm the MeV photons beamed into an angle 1/Tm- The GeV photons can reach us only if 
they are directed in the line of sight, which is along the x-axis. 
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Fig. 2. — A schematic diagram of the geometry of the GeV photon trajectory within the 
MeV shell. The GeV photons meets the MeV photon shell (with a width ~ c^m) at 
Rg and it leaves the shell at -Rc.max- 
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where, in view of the two region model, the MeV component has an upper limit to the 
energy: E M>max . Fortunately t{Eq) is insensitive to the exact value of E M;max . 

At Rm The MeV photons are beamed with an angular width 1/Tm along the radial 
direction outwards. Since R G > Rm the MeV photons travel almost radially outwards. Rel- 
ative to the GeV photons, that move along the x axis, the angular width of the MeV photons 
is of order Rm/{RFm) and it decreases with R. This leads to two angular regimes: Along 
the x axis (for 9 G < R M /(RT M )), the MeV photons have very small angles < Rm/(RTm) 
relative to the GeV photons. This leads to a very small optical depth along the axis. For 
9 G > Rm I (RT m) the angular spread of the MeV photons can be neglected and the typical 
angle between the MeV and the GeV photons is simply R G 9 G /R. 

More generally, the angle between an MeV photon emitted radially outwards at (Rm, 9m) 
and a GeV photon emitted parallel to the x axis at (Rg, 9 g ) i s: 



JJRm sin 9 M ) 2 + (R G sm.9 G ) 2 - 2 (R M sin 9 m )(Rg sin6» G )cos0 

9 = arctan — (5) 

Rg cos 9 g — Rm cos 9m + A 

where A = R — R G is the distance the GeV photon travels before it collides with the MeV 
photon and </> is the angle of a given MeV photon relative to the axis parallel to the GeV 
photon axis. 

A GeV photon, with an (observed) energy E G , and an MeV photon, with an (observed) 
energy Em, can produce a pair if Em > ^M,min> where: 

2(m e c 2 ) 2 4(m e c 2 ) 2 
M ' min ~ (1 + zYE G (l - costf) " (1 + zf9^E G (6) 



The cross section behaves like (e.g. iJauch &: Rohlichlll980[ ): 



^t(1-/? 2 ) 



(3 - /3 4 ) In - 2/3(2 - (3 2 ) 

1-/3 



(7) 



where 



7 = 1/y 1 — /3 2 and a? — 6.65 x 10 25 cm 2 is the Thompson cross section. 
We can estimate now the overall optical depth for the GeV photon: 

t{9 g ,E g )= dR d9 M d4> d E M —^—a{l-cos9) (9) 



Ra Jo Jo JE MtXaia dE M d9 M t 
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The effective optical depth for an observed photon is averaged over all angles: 

exp\-T{E G )\ = — 7, , (10) 



where V = l/(r G (l — /?buik,G cos^g)) is the Doppler factor, /3buik,G — a/1 — l/Tg * s the bulk 
velocity and /?g is the photon index of the GeV emission. It is interesting to note that Eq. 
OH]) is independent of T G . The overall dependence of r on T G arises from Eq. fllOp where 
determines (via T>) the effective width of the integration. 



3. A Simplified model 

We can simplify the above model and obtain an almost analytic formula by making a 
few approximations. We show later that the full numerical results indeed agree with these 
formulae. First, we approximate the MeV spectrum in the relevant energy range using a 
single power law: N MeV (E M ) oc E^ M . Second, we approximate the cross section (Eq. (J7])) 
as a — o"t/(37 2 ) (This form allows us for an analytic integration. The numerical factor is 
chosen by comparison of the analytic approximate results with the full numerical solution.). 
Third, we divide the analysis to two regimes: For 6 G > R M /(RT M ) all the MeV photons 
move radially and the collision angle is simply 9 ~ 9 M . For 9 G < R M /(RT M ) the scatter in 
the directions of the MeV photons is important and we approximate the collision angle as 
9 w R M /{RY M ). 

Define x = R/Rg and x max = Rm^/RG = 1/[1 — 2A M /(9 G R G )]. For large angles the 
collision angle between the MeV and GeV photons is 9 ~ R G 9 G /R = 9 G /x. The minimal 
energy of the MeV photon for pair production is: 



J M,a\m 



4{m e c 2 ) 2 , 



x z . 11 



(l + z)»EaP a 
The number density of the MeV photons: 

77i /i , „\2a„-2r rp a M rp a M~ 2 n 2a M 

N(E M ) * n , MeV (R)(^r a - = i + \l 2( f G ± x-*"»+*>v-°u, (12) 

Ep Att R G c{2iTTn e c ) M 
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where y = Em/Em,^ = 7- Collecting the above expressions: 

dE M -N(E M )(l-cos6)a T /(3f) 



T{V G ) 



dR / 

Ra Je m ^{b,r) 



v l + Zj Lim€q e P CT T U G 



2a M + 1 



6 ■ 4 a Af irR G m e c 3 a, 



1 _ T .-( 2 «M+ 1 ) 
J- J- max 



-(1 + z) 2 e max e G G 



X 1 1 



(13) 



where e G = E G /m e c 2 and e max = E M>max /m e c 2 



We use .Rg = 7 ?2r^cTG/(l + z) to express i?G in terms of T G and r<3, where T G is the 
time when the GeV photon is observed. Generally rj < 1 with 77 = 1 for an external shock 
and 77 < 1(~ 0.01) for an internal shock. As we see later the result is not very sensitive to 
77. Realizing that the GeV photons are mainly coming from 9 G ~ 1 /T G , we can invert now 
this equation to obtain a rough estimate of the minimal Lorentz factor. The solution will be 
consistent if 9 G < Rm/(Rg^m) namely if T G > R g Tm/Rm- To do so we assume that the 
factor in the square brackets of Eq. ( fl3l) is order of unity. We obtain: 



G.min 



(l + z) 2a M-*L M e G 



6 p <J T 



12 • A a M7rm e c 4 7]T G a, 



34 



30 x 10 4q m \ 2-M+ 2 



a 3A Za M 



a,,— 1 a 



^M,51 € GA e P,0 V J G,2 



(14) 



where the notation Q x = Q/10 x is used and ep = Ep/m e c 2 . For a M = 2, the overall 
coefficient is 34. Note the weak (with a power —l/(2a M + 2) ~ —1/6 to —1/2) dependence 
of r Gimin on 77. 

For 6 G < R M /RT M the typical collision angle of the MeV photons is 9 < Rm/(RTm)- 
The optical depth is similar to Eq. f fl3|) . with 6> ~ R M /(RT M ) replacing #g and x max — >■ 00: 



t(#g < i? M /# r 



r(0 G = 0) 
1 



(1 + z y<x M -*L M e 



G 



X 



2« A/ + 1 



6 ■ 4 a M7ri4 aM+1 m e c 3 a M rJ M 



1/1 i \2 / 

-(I + 2) e max e G (— — — . 

4 -Kg-L M 



(15) 



Interestingly, this formula depends on Ym and it is independent of T G . For typical values, 
2 = 1, a M = 2, £ max = ^ = lOGeV, £ P = lMeV, L A/ = 10 51 erg, R M = 10 15 cm, 



R G = 10 16 cm, we obtain T 



Af.min 



8, which provides only a very week constraint on the bulk 
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Lorentz factor of the MeV region. If we set Rg = Rm } this expression reduces to the one 
zone case, where the relevant collisions occur at 6 = 1 /T m in the observer's frame. Taking 
Rm = Rg — 2 T 2 M cST, where ST = 0. 1s is the typical duration of 7-ray pulse, T M > 200, 
consistent with iLithwick fc Saril (120011 ) . 



4. Results 

Fig. [3] depicts the dependence of r on 9g for a set of typical parameters: z — 1, 
St M = 50s and A M = 7.5 x 10 n cm, E G = E max = lOGeV, E P = lMeV, L M = 10 51 erg, 
and T G = 100s. (3g is always chosen to be 2. f ~ 1 for these parameters (with a M = 1, 
Pm = 2) for rG im i n = 40. Note that numerical experiments reveal that to obtain the average 
effective optical depth using Eq. (fTUl) we have to integrate up to ~ 4;/Tg- The results in 
this figure depict both the full (Eq. Q)) and the simplified (Eq. (TT3T) ) calculations. The 
later are depicted by the two thicker lines which are almost superposed on the corresponding 
numerical results, showing consistency of the analytical solution with the full numerical one. 
As the angle increases, the small angle approximation (sin^ ~ 6) breaks down and the 
approximate solution slightly diverges from the numerical one. However, the deviation is 
small and it usually takes place in a regime that is not critical for the overall optical depth. 

One can clearly see different segments of power law dependence of the optical depth on 
9 G - For a given GeV photon 6q determines the colliding angle and hence E M) min- As most 
photons are in the lowest energy and the cross section is largest near EM,min, the optical 
depth is dominated by the low energy photons at i?M,min and the photon index can be 
simply taken the one at -EA/,min- For small angle, E M ,mm > E p (see Eq. and the effective 
spectral slope is (3m- For large angles E M min < E p and the effective spectral slope is a M . 
Correspondingly, lines with the same (3m coincide at small angles while lines with the same 
a M conicide at large angles. These results suggest that the single power law approximation 
for the spectrum is useful. E Mrain determines the relevant spectral, a M or (3 M . 

A second transition takes place when x max approaches unity, namely the width of the 
interaction region is small compared to Rg- According to Eq. ( !T3|) . t(9g) oc 9 g m for 
i max >> 1 and t{6g) oc 6 g m for x max ~ 1 + 2Am / '(9 g Rg) —> 1- The transition between 
the two takes place as 6 G « V / 2A M /R G ~ ^2cT 90 / (2T 2 G cT G ) ~ ^Tw/(T G T G ). For some 
parameters the two transitions may coincide to one and we have chosen Lm so that the two 
transitions are clearly seen. 

Fig. H] depicts r^mm as a function of the MeV luminosity Lm- As one can expect (see 
Eq. P3J) r^mm increases with Lm- The optical depth for Tq = Is (thin solid line, 77 = 1) is 
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Fig. 3. — The optical depth for a photon with Eq = lOGeV as a function of Bq. Different 
lines correspond to different spectral indices. The parameters of the lines are encoded using 
a combination of colors and line types. a M = 0,1,2 are represented as blue, black and red 
respectively. (3m = 1, 2, 3 are represented as dashed, solid and dotted lines respectively. The 
two thick lines are the analytical solution from Eq. ( JTB"]) for a M = (3m — 1 ; 2. 
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much larger than the other ones for which Tq = 100s, again in agreement with Eq. (1141) . If 
a M = Pm, Tcmin increases with L M as a single power law. This is consistent with Eq. (fT4j) . 
For a M 7^ (3m, the relationship breaks into two power law segments that are dominated 
by different parts of Band spectrum of the MeV photons. As the luminosity increases the 
Lorentz factor increases, and the dominant contribution to the optical depth arises from 
6q 1/Tg, with a lower effective 9g, -EA/,min increases and hence for large values of Lm the 
opacity is dominated by the high energy spectral slope, (3 M . Conversely, for low values of 
L M the opacity is dominated by a M . The transition takes place at E M)Ulin cs E p . 

The two zone limits, T^mm, for four Fermi bursts are shown in table [1] together with 
the single zone limits, r m ; n , and some parameters of the bursts. While the single zone 
limits, r min , are of order 1000 and even larger, the two zone limits are around 200 (400), 
for r) — 1 (0.01). It should be stressed that in two of these bursts, GRB 080916c and GRB 
090510, the highest energy GeV photon used to determine the single zone limit is coincident 
with a large MeV flux, while an 11 GeV photon is contemporaneous in GRB 090902b with 
an MeV spike. Still in both GRB 090902b and GRB 090926a GeV photons are observed 
after the end of the prompt MeV (Tg > T 90 ). Thus, it is not clear whether the single zone 
of the two zone limit should be used. 

If the GeV photons are from an external shock, another direct con straint on Tg arise s 
from the dynamics of external shock: Tg — 240EjJ^ 2 rin (1+ z) 3 ^ 8 ^ 3 ^ ' ( Sari Sz Piran 1999 ). 



With observed GeV peak emission time Tq p of ~ 10s (jGhiselrini et al.ll2010f ). E& ~ E 1 
10 55 erg and a circum-burst density n ~ 1cm -3 , we obtain Tq ~ 740. This value is larger 
than the lower limits obtained from the two- zone compactness estimate. However, it is quite 
uncertain, in view of the uncertainty in the determination of T^p- 



Table 1: Limits on Fermi LAT busts 



GRB 


z 


T 90 (s) 




/3m 


-Ep(MeV) 


L iso (erg/s) 


E G (GeV) 


T G 


1 G,min 


Tmin 


ref 


080916C 


4.35 


66 


1.02 


2.21 


1.17 


7 x 10 53 


13.22 


40 


193(414) 


880 


1 


090510 


0.903 


0.5 


0.48 


3.09 


5.1 


4.6 x 10 53 


3.4 


0.5 


150 (277) 


1200 


2 


090902B 


1.822 


30 


0.61 


3.87 


0.8 


3.2 x 10 53 


33.4 


82 


120 (218) 


1000 


3,4 


090926A 


2.1062 


20 


0.693 


2.34 


0.27 


4.2 x 10 53 


19.6 


26 


150 (318) 


1200 


5 



Shown are the low energy spectral parameters a^i, Am and Ep as well as the luminosity, Li SO and the energy, Eg 
and time Tg of the highest energy GeV photon as well as the two zone, TG,mim and the single zone, T m i n , limits. 
(1) Abdo et al. f 2009a') : (2) Ackermann et al. C2010) : (3) Abdo et al. f 2009b) : (4) de Palma et al. (2009') : 
(5) Swenson et al. C20101 

* in bold face is the value for 77 = 1 and brackets is the value for 7/ = 0.01 
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3 




logioL M (erg/s) 

Fig. 4. — The minimal Lorentz factor r^ mm as a function of the MeV luminosity, Lm- 
Different curves correspond to different photon indices. The parameters of the lines are 
encoded using a combination of colors and line types. a M = 0,1,2 are represented as blue, 
black and red respectively. (3m = 1,2,3 are represented as dashed, solid and dotted lines 
respectively. The thick grey line has T 90 = T G = Is, corresponding to a short GRB. 
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Conclusions 



Following various indications that the high energy (GeV) emission in GRBs is produced 
in a different region than the lower energy (MeV) emission we derived here revised compact- 
ness limits on the Lorentz factor of GRB outflows within a two-zone model. We considered 
the "natural" model in which the GeV emission is produced at larger radii than the MeV 
emission. This would arise, for example, if the MeV emission is produced by internal shocks 
and the GeV emission by the afterglow, as has been sugg ested recently by several authors 
(IKumar fc Barniol Duranll2009[ I201CH ; iGhisellini et al.ll2010l ) . We calculated the optical depth 
for pair production by a GeV photon passing through the MeV photons shell. Our results 
reduce to the one-zone model when the emission region of the MeV and GeV coincide. 

Collisions between the GeV and MeV photons occur in the two-zone model at larger 
radii than the prompt emission radius, the density of the MeV emission is smaller (than in 
the prompt emission regime) and the MeV photons are more collimated along the line of 
sight. Consequently, the optical depth is smaller compared to the one- zone case and the com- 
pactness constraint on the Lorentz factor becomes weaker. The new constraint that we find 
is only for the Lorentz factor of the GeV region. The constraint on the Lorentz factor of the 
MeV emitting region, arising from the optical depth of the GeV photon, is rather weak. The 
weak limit does not contradi ct to the neutrino driven ietsdAlov. Janka fc Mullerl I2000T) nor 
to th e magnetic driven jets tlKomissarov et al.l 120091 ; iTchekhovskoy. McKinney fc Narayan 
2009h . 



For a canonical set of parameters, like z — 1, 8t 



M 



50s, E, 



G 



lOGeV, 



E P = lMeV, L M = 10 erg, and T G = 100s, the constraint on the GeV region is T G > 100. 
When we apply the two zone constraint to four Fermi bursts, we find minimal Lorentz factors 
of about 200-400, about one fifth to one half (depending on if) of the the one- zone limit which 
is in the order of 1000. We conclude that one should proceed with care when applying the 
one-zone limits to the Fermi data and unless we can verify that the GeV emission is indeed 
produced in the same region as the lower energy prompt emission we should consider the 
more relaxed two zone limits. 
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